On the effect of weak disorder on the density of states in graphene 
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The effect of weak potential and bond disorder on the density of states of graphene is studied. By 
comparing the self-consistent non-crossing approximation on the honeycomb lattice with perturba- 
tion theory on the Dirac fermions, we conclude, that the linear density of states of pure graphene 
changes to a non-universal power-law, whose exponent depends on the strength of disorder like 
l-'ig/\/^iTt^ , with g the variance of the Gaussian disorder, t the hopping integral. This can result 
in a significant suppression of the exponent of the density of states in the weak-disorder limit. We 
argue, that even a non-linear density of states can result in a conductivity being proportional to the 
number of charge carriers, in accordance with experimental findings. 

PACS numbers: 81.05.Uw,71.10.-w,72.15.-v 

I. INTRODUCTION 

Graphene is a single sheet of carbon atoms with a honeycomb lattice, exhibiting interesting transport 
propertiesi^'^'^'^. These are ultimately connected to the low-energy quasiparticles of graphene, i.e. two-dimensional 
Dirac fermions. Its conductivity depends linearly on the carrier density, and reaches a universal value in the limit of 
vanishing carrier densitj»ii'^. The former has been explained by the presence of charged impurities, while the latter does 
not allow a charged disorderi^. Moreover, in the presence of a magnetic field, the half- integer quantum Hall-effect is 
explained in terms of the unusual Landau quantization and by the existence of zero energy Landau level^ii^. 

The density of states in pure graphene is linear around the particle- hole symmetric filling (called the Dirac point), 
and vanishes at the Dirac point. This is a common feature both in the lattice description and in the continuum. 
In addition, the lattice model also shows a logarithmic singularity at the hopping energy, which is absent in the 
continuum or Dirac description. 

When disorder is present, the emerging picture is blurred. Field-theoretical approaches to related models (quasi- 
particles in a d-wave superconductor) predict a power-law vanishing with non-universal^ or universal^ exponent or a 
diverging^ density of states, depending on the type and strength of disorder. Away from the Dirac point a power-law 
with positive non-universal exponent is also supported by numerical diagonalization of finite size systemsiSiii. At 
and near the Dirac point the behavior of the density of states is less clear. Some approaches favour a finite density of 
states at the Dirac point^^, whereas others predict a vanishing DOS^i^ or an infinite DOSi^. There is some agreement 
that away from the Dirac point and for weak disorder the DOS behaves like a power law with positive exponent 

p{E)^p„\Er (7>0). (1) 

The purpose of the present paper is to investigate, how a non-universal (therefore disorder dependent) power-law 
exponent (found numer ically in Refs. [T^IH 

can emerge for weak disorder (compared to the bandwidth), and what 
its physical consequences are. We determine the exponent based on the comparison of the self-consistent non-crossing 
approximation on the honeycomb lattice and of the perturbative treatment of the Dirac Hamiltonian. The exponent 
decreases linearly with disorder. Then, using this generally non-linear density of states, we evaluate the conductivity 
away from the Dirac point by using the Einstein relation. We show, that based on the specific form of the diffusion 
coefficient, this can result in a conductivity, depending linearly on the carrier concentration^, and in a mobility, 
decreasing with increasing disorder. These are in accord with recent experiment on K adsorbed graphenei^. By 
varying the K doping time, the impurity strength was controlled. The conductivity away from the Dirac point still 
depends linearly on the charge carrier concentration, but its slope, the mobility decreases steadily with doping time. 
Our results apply to other systems with Dirac fermions such as the organic conductor— a-(BEDT-TTF)2l3. 
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II. HONEYCOMB DISPERSION 

We start with the Hamihonian describing quasiparticles on the honeycomb lattice, given by^^^: 

Ho ^ hiai + h2<T2, (2) 
where ct^-'s are the Pauli matrices, representing the two sublattices. Here, 

3 3 
/ii = — < cos(ajk), ft,2 = — i sin(ajk), (3) 

i=i i=i 

with ai = a(— •\/3/2, 1/2), a2 = a(0,— 1) and = a(\/3/2, 1/2) pointing towards nearest neighbours on the hon- 
eycomb lattice, a the lattice constant, t the hopping integral. The resulting honeycomb dispersion is given by 
±y7i|+7i|, which vanishes at six points in the Brillouin zone. To take scattering into account, we consider the 
mutual coexistence of both Gaussian potential (on-site) disorder (with matrix element Vg^ri satisfying (Vo,r) = and 
variance (K),rVo,r')5o = Qo^rr') and bond disorder in only one direction (in addition to the uniform hopping with 
matrix element Vb^n satisfying (V^.^) = and variance {Vb^rVb.r') = gtSrr'), which is thought to describe reliably the 
more complicated case of disorder on all bonds^^. In graphene, ripples can represent the main source of disorder, and 
are approximated by random nearest-neighbour hopping rates, while potential disorder might only be relevant close 
to the Dirac point^. The corresponding term in the Hamiltonian is 

V ^Vo,r<JO + Vb^r<Jl, (4) 

which results in H = Hq + V. 



FIG. 1: (Color online) A small fragment of the honeycomb lattice is shown. The thick red lines denote the uni-directional bond 
disorder, on-site disorder acts on the lattice points. 

Without magnetic field, the self-energy for the Green's function, which takes all non-crossing diagrams to every 
order into account (non-crossing approximation, NCA). can be found self-consistently from^S 

1 - (go + gb)G^ ^ 

[igo + 9b)-i9o-9brG^]G 
where iujn is the fermionic Matsubara frequency and 

G = Go[it^„ - S(iw„)]. (6) 
Here, Go is the unperturbed local Green's function on the honeycomb lattice given by 

Q / ^ _ f ZCpk 

"^""^ ~ (27r)2 7 z2-t2[4cos(V3fc:,/2)cos(3fcy/2)-h2cos(V3/c^)-h3]' 

where Ac ~ 3\/3a^ /2 is the area of the unit cell, and the integral runs over the hexagonal Brillouin zone with corners 
given by the condition + /ij = 0. This can further be brought to a closed form using the results of Ref. [2l|. On the 
other hand, in the continuum representation, using the Dirac Hamiltonian, the above Green's function simplifies to 

^ , 2Acao f (fk AcZOQ , / A2\ 
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and V = 3ta/2. The cutoff A can be found by requiring the number of states in the Brillouin zone to be preserved in the 

Dirac case as well. This leads to A = \/ iry^t. Another possible choice relies on the comparison of the low frequency 
parts of the Green's function in the lattice and in the continuum limit, which reveals the presence of ln(tj/3i) terms. 
This leads to A = 3t, which coincides with the real bandwidth on the lattice. We are going to use this form in the 
following. The difference of the variances becomes important when calculating the 2nd order correction (in variance) 
to the self-energy. It is clear from Eq. ([5]), that the same self-energy is found for pure potential or unidirectional bond 
disorder. The effect of their coexistence is the strongest, when they possess the same variance. From this, the density 
of states follows as 

p{lj) = --ImGicj + ie) (9) 

TT 

with e ^ 0+. Without disorder, wc have the linear density of states p{uj <C t) = Ac\uj\/27tv'^. At zero frequency, in 
the limit of weak disorder, the sclf-cncrgy is obtained as 

S](0)^-»Aexp('- ), (10) 

which translates into a residual density of states as 

A / TTV^ \ , , 

pO = ^ ^exp . 11 

From this expression, weak disorder is defined by the condition go + g^ <C ■ The exponential term indicates the 
highly non-perturbative nature of density of states at the Dirac point: all orders of perturbation expansion vanish 
identically at w = 0. 

From Eq. ([5]) it is evident, that the interference of the mutual coexistence of both on-site and bond disorder should 
be the most pronounced when go — gt- The frequency dependence of the density of states on the honeycomb lattice 
can be obtained by the numerical solution of the self-consistency equation, Eq. ([5|), and is shown in Fig. [2l For small 
frequency and disorder, there is hardly any difference between pure on-site or unidirectional bond disorder and their 
coexistence. However, at higher energies and disorder strength, they start to deviate from each other. At lo = t, 
the weak logarithmic divergence is washed out with increasing disorder strength. Such features are absent from the 
Dirac description, which concentrates on the low energy excitations. Interestingly, for weak disorder, the residual 
DOS remains suppressed as suggested by Eq. (jlip . but the initial slope in frequency changes. In order to determine, 
whether the exponent or its coefficient or both change with disorder, we perform a perturbation expansion in disorder 
strength using the Dirac Hamiltonian to quantify the resulting density of states, and compare it to the numerical 
solution of the self-consistent non-crossing approximation using the honeycomb dispersion. 



III. POWER-LAW EXPONENT 



The expansion of the one-particle Green's function in disorder at z = E + ie leads to 

G{z) = Go + GoVGo + GnVGoVGo + 



(12) 



where V = K/r^o -I- Vb,rcri describes both Gaussian potential and unidirectional bond disorder, G ~ {z — Hq — V) ^ 
and Gq = {z — Hq)^^, and Hg = v^k^ui + ky(T2) is the Dirac Hamiltonian. After averaging over disorder, we get 



(Grr) — Go.rr +.9 ^ Gp-rr' Gp-r'r' Gp-r'r 



(13) 



with g = go + gb- The Green's function Go is translational invariant and reads from Eq. ([8]) at real frequencies as 



Go;?- 



Ac\E\ 
2ttv^ 



(14) 



This implies 



(Grr) = Go;rr + giGl)rrGo;rr + o{g'^) = Go;rr[l + 5(G'o)rr] + o(g^) 



(15) 
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{go + 9b)/t' 



FIG. 2: (Color online) The density of states is shown in the left panel for pure on-site [gb = 0) or unidirectional bond [go = 0) 
disorder (solid line) for {go + gb)/t^ ~ 0.1, 0.4, 0.8 and 1.6 with decreasing DOS ak, uj — t. The red dashed line represents 
the coexisting bond and unidirectional bond disorder with go = gb- The black dashed-dotted line denotes the free case with a 
linear density of states at low energies, exhibiting a logarithmic divergence at tj = t in the pure limit. The right panel shows 
the residual density of states for on-site or unidirectional bond disorder (blue solid line) and their coexistence with go — gb (red 
dashed line). The black dashed-dotted line denotes the approximate expression, Eq. (|ll|l . for weak disorder. For go + gb < OAt^ , 
the residual density of states is negligible. 
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FIG. 3: (Color online) The low energy density of states is shown in the left panel for {go -f gb)/^^ ~ 0.004, 0.1, 0.2, 0.3, 0.4 and 
0.5 from bottom to top, for pure on-site [gb = 0) or unidirectional bond {go — O)disorder at positive energies (blue solid line), 
and for their coexistence at negative energies at go = gb (red dashed line) . The black dashed-dotted line denotes the power- law 
fit as p{uj) = po + 2pi{ui/t)'^ . The green vertical dotted line separates the two parts. The right panel visualizes the exponents as 
a function of the variance of the disorder, go+gb, for on-site or bond disorder (blue solid line) and their coexistence {go = gb)- 
The black dashed-dotted line denotes the result of perturbation theory: 7 = 1 — ^{go + gb) / ^/ii^t^ - As is seen, the agreement 
is excellent in the limit of weak disorder. Note, that g denotes the variance of the disorder. 



Moreover, we have with Eq. ([14]) 



(27ri;2)2 7 [z + v{k^ai + kya2)Y 



dz 



ln(l - 



2A2 



^2 _ ^2 



(To 









— m 







(16) 
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for \ ^ \E\ e. Therefore, we obtain 

(Grr) = Go:rr 

From this, the density of states follows as 



p{E) = lm{Gr 

n 







igAc' 






2z;2 _ 



0(5^)- 



AcE 



(17) 



(18) 



If we further assume that (I) the density of states as a function of E satisfies a power law (inspired by Refs. lol fiolflll ) 
and (II) the disorder strength g is small, we can formally consider Eq. (jlSp as the lowest order expansion in disorder, 
and sum it up to a scaling form as 



p{E) 



AcX 



l-(2gAc/7TO^) 



(19) 



Hence, this suggests that the linear density of states of pure graphene changes into a non-universal power-law depending 
on the strength of the disorder as 7 ~ 1 — {Ag/n^/St'^). Note, that the exponent does not depend on the ambiguous 
cutoff A. We mention that Eq. (jl8|l might suggest other closed forms than Eq. p9|) . By using the renormalization 
group procedure to select the most divergent diagrams at a given order g (similarly to parquet summation in the Kondo 
problem), one can sum it up as a geometrical seriesi^iS^. However, the resulting expression contains a singularity 
around |S(0)| (Eq. (fTO|) . playing the role of the Kondo temperature here), and is valid at high energies compared 
to_[S(0)| as Eq. (fT8|) . To avoid such problems, we use a different scaling function, suggested by the results of Refs. 

We compare this expression to the numerical solution of the self-consistent non-crossing approximation in Fig. [S] 
To extract the exponent, we fit the data with p{E) = pQ + 2pi{\E\/t)"' , and extract po.i and the exponent 7. As can be 
seen in the left panel, the power-law fits are excellent in an extended frequency window up to t/A. This suggests that 
this effect should also be observable experimentally as well. The obtained value of po is negligibly small, as follows 
from Eq. pT|) . From the fits, we deduce the exponent and its coefficient, which is shown in the right panel. It agrees 
well with the result of perturbation theory, Eq. (|19p in the limit of weak disorder. The suppression of the exponent is 
significant, and can be as big as 30-35% around {go + gbj/t^ ^ 0.4. Similar phenomenon has been observed for Dirac 
fcrmions on a square lattice in the presence of random hoppin g^^i^^ , where disordered systems were studied by exact 
diagonalization. The decreasing exponent with disorder agrees with our results. 



IV. CONDUCTIVITY FOR NON-LINEAR DENSITY OF STATES 

Now we turn to the discussion of the conductivity in graphene. A possible starting point is the Einstein relation^^, 
which states for the conductivity 

CT = e^pD, (20) 

where p is the density of states and D is the diffusion coefficient, both at the Fermi energy Ep- Assuming a general 
power-law density of states, as found above, we have p{E) — pi{E/X)'^ . In the weak-disorder limit and away from the 
Dirac point E' = 0, we can safely neglect any tiny residual value. Moreover, the diffusion coefficient in this case is of 
the form D = DiE, which is validated from the Boltzmann approach in the presence of charged impurities^. On the 
other hand, at the Dirac point there is a exponentially small density of states and a finite non-zero diffusion coefficient 
D oc g/piii the presence of uncorrelated bond disorder— such that the conductivity is of order 1, in units of e^/h . In 
the following, however, we will concentrate on the regime away from the Dirac point. Putting these results together, 
we find 

^7+1 

cT = e'piDi^. (21) 

This can be simplified further by noticing that the total number of charge carriers, participating in electric transport, 
can be expressed as 

^ p{E)dE ^ p,j,:^f^. (22) 
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By inserting this back to Eq. ([2T|) . wc can read off the conductivity as 

a ^ e^Diij +l)n. (23) 

From this we can draw several conclusions. First, it predicts that away from the Dirac point, where our approach 
predicts a general power-law density of states, the conductivity varies linearly with the density of charge carriers, in 
agreement with experiments^. Second, the mobility of the carriers, which is the coefficient of the n linear term in the 
conductivity, behaves as 



where we used our approximate expression for the exponent in the density of states, Eq. (|19[) . This means, that 
with increasing disorder (g), the mobility decreases steadily, in agreement with recent experiments on K adsorbed 
graphenei^. There, the graphene sample was doped by K, representing a source of charged impurities. Nevertheless, 
these centers also distort the local electronic environment, and act as bond and potential disorder as well. The 
observed conductivity varied linearly with the carrier concentration n, similarly to Eq. ([23)) . Moreover, the mobility 
(the slope of the n linear term) decreased steadily with the doping time (and hence the impurity concentration), 
which, in our picture, corresponds to a reduction of the exponent 7 as well as the mobility, Eq. (|24p . 

To study the properties close to the Dirac point we have to go beyond the perturbative regime. Then we realize 
that the density of states does not vanish at E = 0. As an approximation we add a small contribution near the Dirac 
point 



piE) ^ poS„{E) + pi 

in form of a soft Dirac Delta function 



IT E^ + ■n'^ 



77- 

This implies a particle density n which does not vanish at the Dirac point: 



niEp) = / p{E)dE « po + T^^^^F • (25) 



/o ' ' ' (7 + l)AT 

Moreover, the diffusion coefficient does neither diverge nor vanish at the Dirac poin t^^i^^ such that we can assume 

D{E) = DoSr,{E)+DiE . 

From the Einstein relation we get the conductivity which provides an interpolation between a behavior linear in n 
away from the Dirac point and a minimal conductivity at the Dirac point: 



(26) 



DapaSf^iEp) for Ep ^ 

DipiEj^^/X'' - (1 + 7)71 for £;f > 



This, together with Eq. (|25p . implies for Ep ^ the same behavior as in Eq. ((23)) with the mobility of Eq. ([24 
The value of the minimal conductivity can be adjusted by choosing the parameter 77 properly. Therefore Eq. 
provides us with a qualitative understanding of the conductivity in graphene for arbitrary carrier density. 



V. CONCLUSIONS 



Wc have studied the effect of weak on-site and bond disorder on the density of states and conductivity of graphene. 
By using the honeycomb dispersion, we determine the self-energy due to disorder in the self-consistent non-crossing 
approximation. The density of states at the Dirac point is filled in for arbitrarily weak disorder. Wc investigate 
the possibility of observing non- linear density of states away from the Dirac point, motivated by numerical studies 
on disordered Dirac fermionic systems. By comparing the results of non-crossing approximation on the honeycomb 
lattice to perturbation theory in the Dirac case, we conclude, that a disorder dependent exponent can account for 
the evaluated density of states. The exponent decreases linearly with the variance for weak impurities. Then, by 
using the obtained power-law DOS, we evaluate the conductivity away from the Dirac point through the Einstein 
relation. Wc find, that this causes the conductivity to depend linearly on the carries concentration by assuming that 
the diffusion coefficient is linear in energy^, and the mobility decreases steadily with increasing disorder. These can 
also be relevant for other systems with Dirac fermionsi^. 
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